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1. Introduction 

The simplest and perliaps tlie only possible microscopic realization of the two- 
dimensional quantum gravity is provided by the ensemble of planar graphs |jl|]. In its 
turn, the ensemble of planar graphs is generated by the (divergent) quasiclassical expan- 
sion of N X N matrix integrals, with the parameter playing the role of topological 
coupling constant. The formulation of 2D quantum gravity in terms of random matrix vari- 
ables opened the possibility to apply powerful nonperturbative techniques as the method 
of orthogonal polynomials [0 and led to the discovery of unexpected integrable structures 
associated with its different scaling regimes. Originally, a structure related to the KdV 
hierarchy of soliton equations was found in the continuum limit of pure gravity The 
KdV hierarchy can be obtained as a 2-reduction of the KP hierarchy, which defines the 
integrable structure of the microscopic theory. Consequently, the partition functions of the 
A series of models of matter coupled to 2D gravity were identified as r-functions of higher 
reductions of the KP integrable hierarchy Similar statement concerning the D series 
was made in 0. 

The description of interacting gravity and matter based on the higher reductions of 
the KP hierarchy is well suited for studying the spectroscopy and the correlation func- 
tions of local scaling operators. On the other hand, it becomes less efficient for evaluating 
macroscopic loop correlators, topology changing interactions and other infrared phenom- 
ena. Indeed, at distances much larger than the correlation length of the matter fields, the 
relevant integrable structure is this of pure gravity. Therefore it make sense to look for 
a complementary description in which the matter is considered as a perturbation of pure 
gravity. The simplest realization of this idea is given by the Kazakov's multicritical points 
of the one- matrix model P] , which describe a series of nonunitary matter fields coupled to 
gravity. 

A systematic construction of the matter fields as deformations of pure gravity will be 
the subject of this talk, which is partially based on the work 0. Our main task will be 
to exhibit in any theory with matter fields the integrable structure inherited from pure 
gravity. Namely, we are going to derive bilinear functional equations that can be viewed 
as a deformation of the Hirota equations for the partition function of pure gravity. We 
are going to consider in details only the microscopic realization of the theory, but our 
construction survives without substantial changes in the continuum limit. 

We start with the derivation of the Virasoro constraints and the Hirota bilinear equa- 
tions for the matrix integral generating the ensemble of planar graphs. The Virasoro con- 
straints follow from the loop equations, which relate boson-like quantities (traces), while 
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the Hirota equations follow from orthogonality relations involving fermion-like quantities 
(determinants). Then we show how both types of equations generalize for the matrix in- 
tegrals describing theories of gravity with matter fields. For this purpose we will exploit 
the microscopic construction of 2D quantum gravity given by the random-lattice versions 
of the s/(2)-related statistical models: the Ising Q and the 0(n) models, the SOS and 
RSOS models and their ADE and AbiS generalization. H-QB. Eaclr of tlrese models 
can be reformulated as theory of one or several random matrices with interaction of the 
form tr ln(Ma 01 + 1® Mi). Such an interaction can be introduced by means of a differen- 
tial operator of second order H acting on the coupling constants; the partition function of 
the model is obtained by acting with the operator on the partition function of one or 
several decoupled one-matrix models. The operator defines a canonical transformation 
of the symplectic structure associated with the space of coupling constants. 

The Virasoro constraints = for each of the one-matrix integrals transform to 
linear differential constraints e^L'^e~^ = for the interacting theory. These constraints 
are equivalent to the loop equations, which have been already derived by other means. 
The new point is that the Hirota bilinear equations that hold for each of the one-matrix 
integrals induce, by simply replacing the vertex operators as V±(z)" — > e^V^(2)e~"^, 
bilinear functional equations for the interacting theory. In this way the integrable struc- 
ture associated with pure gravity is deformed but not destroyed by the matter fields. In 
particular, the deformed Virasoro and vertex operators satisfy the same algebra as the 
bare ones. 

In the continuum limit, the loop equations and the bilinear equations are obtained as 
deformations of the Virasoro constraints and the Hirota equations in the KdV hierarchy. 
Their perturbative solution is given by the loop-space Feynman rules obtained in . 



2. The one-matrix model 

The partition function of the ensemble of all two-dimensional random lattices is given 
by the hermitian N x N matrix integral 

/oo 
dM exp (^tvJ^tnM'^y (2.1) 

n=0 



^ The g-state Potts |12] needs a more elaborate construction and will not be considered here. 
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The integrand depends on the matrix variable M only through its eigenvalues A^, i = 

1.. ..,N. Retaining, therefore, only the radial part of the integration measure dM ~ 
rii^i '^-^i rii/jl-^i — Xj) we can write the integral ( |2.1|) as the partition function for a 
system of Coulomb particles in a common potential 

ZN[t] = I l[dX, exp (^t„A7) 1[{X, - Xjf. (2.2) 

i=l n=0 i<Cj 

Let us recall that the last factor is the square of the Vandermonde determinant 

An{Xi,...,Xn)= n (Az-A,) = det[(A,)^-i]. (2.3) 

l<i<j<N 

In the following we will denote by (...) N,t the mean value defined by the partition function 
( pj.2| ). For the time being we restrict the integration in A's to a finite interval [A^, Ajj] on 
the real axis, so that the measure 

oo 

dfit{X) = rfAexp (^^t^A'^) (2.4) 

n=0 

is integrable for any choice of the coupling constants tn- The limit Xl —oo^Xr oo 
exists if the exponential vanishes at ±oo. 

2.1. Loop equations and Virasoro constraints 

The loop equations determine the expansion of the integral (|2.2|) and represent an 



infinite set of identities satisfied by the correlation functions of the collective loop variable 

N 

1=1 * 

The derivation goes as follows. From the translational invariance of the integration measure 
dXi we find, neglecting the boundary terms at A = Xl,r, 

( E (I- + 2 5: ^ ^ 5: n*..AT') r^)„ , ^ 0. (2.6, 

i=l n>0 

Using the identity 

V ^ I 2 V ^ ^ - V ^ ^ (27) 
^{z- AO^ ^.z-\X,-Xj ^z-X,z-Xj ^ ■ ^ 
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we find 

N 



1=1 * n>0 



The sum in the second term can be viewed as the result of a contour integration along a 
contour C enclosing the interval [A^, Ajj] but leaving outside the point z: 



Introducing the collective field 



n>0 

we write the last equation as 



^ =0. (2.11) 



1(2 27rz 2 — 2' \ / Af,t 

The insertion of the operator trM"^ = ^? '^^^ t)e realized by taking a partial 
derivative with respect to the coupling The collective field ( |2.10| ) is therefore represented 
by the free field operator acting on the coupling constants 

*(.,4E'-"-1'>(i)^-E^I: 

n=0 n=l 

and the loop equation ( p.ll| ) is equivalent to the linear condition 

f dz' 1 

i — -T{z') ■ Z^[t\ ^ {T{z) . Z^[t\)^ = (2.13) 

imposed on the partition function, where 

T{z) = \ : {d^{z)f : (2.14) 
denotes the energy-momentum tensor for the gaussian field (|2.12| ). Expanding 

: (a*(/))':=5^L^^-""' (2.15) 

we find a set of differential equations (Virasoro constraints) 



L„-ZAr[t] = (n>-l) (2.16) 
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where 

Ln = E7^757^ + E^^^7^ (2-17) 
fro ^^-^ fro ^^-+^ 

satisfy the algebra 

[L^,L„] = (m - n)L^+n- (2-18) 
With respect to the trivial variable to, we have the additional equation 

-^ZN[t] = NZN[t]. (2.19) 

Oto 

An elaborate analysis of the loop equations and explicit expressions of the lowest 
orders of the 1/N expansion of the free energy are given in ref. [H . 



2.2. Orthogonal polynomials and Hirota equations 



It is known [|T5[ that the partition function (|2.1|) is a r-function of the KP hierarchy 
of soliton equations. The global form of this hierarchy is given by the Hirota's bilinear 
equations [0 (for a review on the theory of the r-functions see, for example, fl^.) Below 



we will derive the Hirota equations using the formalism of orthogonal polynomials. 
Let us define for each {N = 0, 1, 2, ...) the polynomial 

P^,t(A) = (det(A-M));v,t 

N 

= (n(^-^^))^, (2-20) 
= A^ - {trM))N,t A^-^ + ... + (-)^(det M))Ar,t. 

It easy to prove that the polynomials ( 2.20| ) are orthogonal with respect to the measure 
(i/Ut(A). Indeed, 

Zn J dHt{XN+l)PN,ti^N+l)^N+l 

.N+i (2.21) 

/ Yl dHt{X^)AN+liXl,...,XN+l) 



i=l 



N+1 



(-)^+^ ''A^AAr+i(Ai, As, Xn+i)- 



The sum in the integrand is the expansion of the determinant 



Ai ■ 
A2 ■ 


■ Af-1 

\N-1 


A^ 
Ai 


N+l ■ 


\N-1 





1 A; 

with respect to its last column. It vanishes for k — 1, N — 1, which proves the statement. 



For n — N, one finds 



Hence 



/ d//t(A)P,v,t(A)A^ = Zn+i/{N + 1). 



(2.22) 



(2.23) 



'A, iN + l)ZN[tY 

The orthogonality relations ( p.23|) can be written in the form of a contour integral, 
namely, 



>N,k 



(2.24) 



where the integration contour C encloses the point z = and the interval [A^, Ajj]. Indeed, 
the residue of each of the n poles is equal to the left hand side of ( p.23| ) multiplied by 

Zk[t]/Zk+i[t]. 

A set of more powerful identities follow from the fact that the polynomial PN,t{z) is 
orthogonal to any polynomial of degree less than and in particular to the polynomials 
Pk,t' [A], /c = 1, 2, A — 1 where t' = {t^, n = 1, 2, ...} is another set of coupling constants. 
Written in the form of contour integrals, these orthogonality relations state, for A' < A, 



e^r=i(*"-*")^7 det(2-M) ) ( — — ^ r 



0. 



(2.25) 



N',t' 



The Hirota equations for the KP hierarchy are obtained from eq. (|2.25|) after express- 
ing the mean values of (det^^) in terms of the vertex operators 



^^±{z) = exp ± ^ exp T In - 7^ T J] ^ 



n=0 



z dtn ^ n dtr 

n=l 



It follows from the definition ( |2.2| ) that 

V±(z) • ZN[t] = e^^r=o*"^"(^det(z-M)±^^^^Z7v[t] 
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(2.26) 



(2.27) 



and eq. (|2.25|) is therefore equivalent to 

j^dz (\+{z)-ZN[t]) (\-{z)-ZN'[t']) =0 

which is one of the forms of the Hirota equation for KP [jl^ . 
After a change of variables 



(A^' < N), 



(2.28) 



t +t' 



Vn 



t -t' 



2 ' " 2 ' 

the Hirota equations (|2.28| ) take its canonical form 

Res^=o z^-^' e^^^i'^y^^^^ e~^^=i^^ Zn[x + y]ZN'[x - y] =0 



(2.29) 



(2.30) 



where N' < N . The differential equations of the KP hierarchy are obtained by expanding 
( p.30|) in y^. For example, for N' = N, the coefficient in front of yf is 



+ 3 



-4- 



d d 



dyf dyi dyi dys 



ZN[t + y]ZN> [t - y] 







y=o 



o2 

and one finds for the "free energy" u[t] = 2^ logZjv the KP equation 



, du „ du d^u 
ot'i dti dt{ J 



0. 



(2.31) 



(2.32) 



The lowest equation in the case N' = N ~ 1 (modified KP) is the so called Miura trans- 
formation relating the functions u = 2df log ZN[t] and v = log(Zjv+i [t]/-Z'Ar[t]): 

u = d2V- dfv - {divf. (2.33) 

One can check directly that these differential equations are satisfied by the asymptotic 
expansion around the gaussian point 



ZN[t] = {-t2) / exp(^Ato-— - + — ^--^ + 



^2.34) 



The Hirota equations themselves are not sufficient to determine uniquely the partition 
functions -Z^Ar[t]. The additional information comes from the "string equations" 



y^nt^ [ dllt{\)\^-^PN,t{\)PN-l+a,t{\) = ^0,^4^ (^ = 0,1) 



(2.35) 



obtained, with integrating by parts, from 

/ d^t{X)PN-l+a,tWdxPN,tW = Sa,o4^ = 0' 1) (2-36) 

and therefore valid only if dfj,{X)/dX vanishes at the endpoints Xl^Xr of the eigenvalue 
interval. 

In terms of vertex operators, the string equations (|2.35| ) read 

Vnt, / dz z^-^ (V_(^) ■ Z^+,[t])(V+(^) ■ ZN[t]) = 6^,oN Z^ltf (2.37) 

with cr = 0, 1. 
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3. The 0(n) model 

The partition function of the 0(n) matrix modeli is defined by the N x N matrix 
integral 

n=0 n+m>l 

(the sum over n and m runs the set of nonnegative integers) and describes the ensemble of 
nonintersecting loops on a random graph. A loop of length n is weighted by a factor nT""^ 
where T is the temperature of the loop gas. For nG [—2,2], the model exhibits critical 
behavior with spectrum of the central charge C = 1 — 6{g — 1)^/(7, 9 = ^ arccos(— n/2); the 
critical behavior for |n| > 2 is this of a branched polymer [0. The matrix integral ( |3.1j ) 
can be again interpreted as the partition function of a Coulomb gas with more complicated 
but still pairwise interaction between particles: 

^o,„,,^, ^ / n^A. exp (f «„A?) ^Z'l^\,^l (3-2) 

The interval of integration [A^, Xji] should be such that < T/2. With the last restriction 
the denominator in the integrand never vanishes. The choice of the integration interval 
does not influence the quasiclassical expansion and hence the geometrical interpretation in 
terms of a gas of loops on the random planar graph. The saddle point spectral density is 
automatically supported by an interval on the half-line [— oo,T/2] 0. The most natural 
choice for the eigenvalue interval is therefore Al — > —oo,Xr T/2, with the conditions 
dfi{X)/dX\-^ = d^{X)/dX\T/2 = 0. 

The 0(n) model reduces to the hermitian matrix model in the limit n— > and/or 
T ^ oo, and can be considered as a deformation of the latter in the following sense. Let 
us define the differential operator 



2 



dtn ^ n + m nl ml dt^ dt^ 

n+m>l 



(3.3) 



acting on the coupling constants. It is easy to see that the partition function ( p.2|) is 
obtained from the partition function of the one-matrix model by acting with the operator 



Z2(")[t]=e"-Z^[t]. (3.4) 



^ We use a Roman letter for the parameter n€ [—2, 2] to avoid confusion with the summation 
index n running the set of natural numbers. 
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This simple observation will be of crucial importance for our further consideration. It 
means that the integrable structure of the one-matrix model survives in some form in the 
0{n) model. The operator ( ^.3| ) defines a canonical transformation 

d d 



E^_;fi_^ + — 1)! d 



(3.5) 



n\ ml dtm 
preserving the symplectic structure. 

3.1. Loop equations and Virasoro constraints 

From the translational invariance of the integration measure dXi we find 

i^jhif ,4^^-^ - E 74^-"^r')„,, - «■ (3.6, 

or, expressing the sum as a contour integral, 

Wizf+(( ^^—Wiz')[nW{T-z') + yntn{zT~^]) =0. (3.7) 
Jr 2m z — z ^-^ I N,t 

The integration contour C_ encloses the interval [Xl^Xr] and leaves outside the interval 
[T — Xn, T — Xl] and the point z. Introducing the collective field 

$(z) = tr log{z - M) - -tr log(T - z - M) + ^nz'^ (3-8) 



n 



we write the last equation as 



2'Ki z — z' \ / N,t 

The field ( |3.8| ) is represented by the linear operator 



^ ,({d^z'))')^^ =0. (3.9) 



. i E - (^) I - E --^f-^- #- (3.10, 



2 ^ \ z / dtn ^ n dt. 

n=0 n=l 

and is related to the "bare" field ( |2.12|) by *(z) = e^^{z)e~^. The loop equation 
is therefore equivalent to the linear differential constraints 

L^-Z^(")[t] =0 (n>-l) (3.11) 

where the operators 

~ >^ d d sr^ , d 

= 2^ BrTTf — + ^ ^^''TTf — 

f— ' Otk Otn-k ^ Ol'n+k 

k-m-i{k + m)\ d d 



+ n$:T 



dt 



are related to the standard Virasoro generators ( p.l7|) by = e^L^e"^ and therefore 
satisfy the same algebra. 
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3.2. Bilinear functional equations 



In quite similar way the Hirota equations ( |2.28| ) provide, due to the relation 
set of bilinear equations for the partition function of the 0(n) model 



dz(V^{z)-Z^^-\t]) (V_iz)-Z)^\->[t 

where y±{z) are the transformed vertex operators 
V±=e^V±{z)e-^ 

/ oo 



.0(n) r,/i 







(3.13) 



exp =1= In 



(T - zr 



n=l 
oo 



(3.14) 



n=l 



2 — n(T — z) 



n 



_d_ 

dt„ 



After being expanded in tn ~ t'^, eq. ( |3.13| ) generates a hierarchy of differential equations, 
each of them involving derivatives with respect to an infinite number of "times" tn- 
The functional relations (|3.13| ) are equivalent to the bilinear relations 



dz 



/ det(^ - M) \ / det(T -z- M)" 
\det{T-z-MY / ^\ det(z - M) 



(3.15) 



(A^' < A^) 



[N',t'] 



where ( ) N,t denotes the average corresponding to the partition function (|3.2| ) . These 
relations can be also proved directly by exchanging the order of the integration in the A's 
and the contour integration in z. The "string equations", needed to determine completely 
the partition function, are obtained from ( p.37|) by replacing \± — > y^±{z). 



4. ADE and AdE models 

The ADE and ADE matrix models give a nonperturbative microscopic realization 
of the rational string theories with C < 1. Each one of these models is associated with 
a rank r classical simply laced Lie algebra (that is, of type Ar, Dr, Eqj^s) or its afine 
extension, and represents a system of r coupled random matrices. The matrices Ma of size 
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Na X Na (a = 1, 2, r) are associated with the nodes of the Dynkin diagram of the simply 
laced Lie algebra, the latter being defined by its adjacency matrix G with elements 

Qab _ 1 1 if the two nodes are the extremities of a link < ab > ^-j 
1 otherwise. 

The partition function Z^[i] depends on r sets of coupling constants t = {t^ \ a = 
l,...,r; n = 1,2,...}. The interaction is of nearest-neighbor type and the measure is a 
product of factors associated with the nodes a and the links < a6 > of the Dynkin diagram 



P r / r oo 



+ E ^^^trM^trM,"" 

2^ ^m + 7im!n! 

a, fa m,n=l 

The target space of the Ar model is an open chain of r points and the one of the Ar-\ 
model is a circle with r + 1 points. In this sense the 0(2) model can be referred to as 
the model of the A series. The continuum limit of the Aj. model is that of a free field 
coupled to gravity and compactified at the radius r in a scale where the self-dual radius is 
r = 2. 

Again, the only nontrivial integration is with respect to the eigenvalues \ai [i — 
1, No) of the matrices M^: 



^ ,n , 71, I I _ I I . • 1 Ao,- — Aoi n 

n<ab> U.i,jiT - Xai - Xbj) 



a=l 1=1 

The domain of integration is assumed to be a compact interval [Al, Xr] with Xr <T/2. 

The partition function ([4.3| ) can be obtained by acting on the product of r one-matrix 
partition functions ZN^[t°'],a — 1, ...,r, with the exponent of the second-order differential 
operator 



namely. 



2 

a,b 



^ dt?. dtk ^ n + m n\ m\ dt"; dtt, 

" 'J n+m>l 



(4.4) 



4[^"]=^"-n^^^[^i- (4-5) 



a=l 
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4.1. Loop equations and Virasoro constraints 

The Dyson-Schwinger equations are written in terms of the r-component collective 
variable 

Wa{z) = Y,^—. (4.6) 

1=1 



They are obtained as in the previous cases, by shifting the variable Xai, and read 
dz' 1 

b 

The loop equations are equivalent to the differential constraints 



Wa{zf+i + =0. (4.7) 

JC- ^T^l Z — Z — ^ ^ — ^ J / N,t 



t^-Z^[t]=0 (n> -1, a = l,...,r), (4.^ 



where the operators 



d d d 

fc=0 ^ ^^n-k '-'^n+k 

.sr^Qabsr- r^-k-m-i ik + m)l d d 
6 fc;: dt^n+k+idtl 



(4.9) 



are related to the "bare" Virasoro generators by = e^L^(2;)e ^ and form r commuting 
Virasoro algebras 

[t^^,i':,]=6aA^-n)t^^+^. (4.10) 
With respect to the trivial variables to, we have the additional equations 

— Z^[t]=N,Z^[t], (a = l,...,r). (4.11) 

4.2. Bilinear functional equations 

Each of the one-matrix partition functions on the right hand side of ( ^.5| ) satisfies the 
Hirota equations (|2.28D . As a consequence, the left hand side satisfies a set of r functional 



equations associated with the nodes a = 1, ...,r of the Dynkin diagram. Introducing the 
transformed vertex operators 

Yl{z) = e^Vl{z)e-'^ (4.12) 
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we find, for eacfi node a = 1, r, 



j( ■ Z'^[t\) (V^(^) ■ Z%[i']) =0 {K< iV„) (4.13) 



where tlie integration contour C_ in ( [4.13| ) encloses tlie interval [Xl, Xr] but leaves outside 
the interval [T — Ai?, T — A^]. 

Inserting ( [4.4| ) in the definition ( 4.12| ) we find the explicit form of the dressed vertex 
operators 



n(^) 

JJ(T-22)-5G"''e^^rLo*"^"exp ( t 
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exp 



n=l 



(4.14) 



n=l 



and write the generalized Hirota equations ( [4.13| ) in terms of correlation functions of de- 
terminants: 



dz 



det(z - Ma] 



c_ Y[,{T - 2z)G-' \ n, det(T - z - M,)^- 



nfcdet(r-z-Mb 
det(;z - Ma) 



(4.15) 



{K < Na). 



[N',t'] 



Here ( )g y denotes the average in the ensemble described by the partition function ([4.3|). 
It is possible to prove eq. ( |4.15|) directly by performing the contour integration. It is 
essential for the proof that the interval [T — Xr,T — Xl] is outside the contour C_. 



5. Continuum limit 

We conclude with a brief description of the continuum limit of the above construction. 
In the scaling limit the field ( |2.12| ) is expanded in the half-integer powers of the (shifted 
and rescaled) variable z 

dHz) = \'£t.z^-Y.'-^ir (.G2 + 1) (5.1) 

s>0 s<0 ^ 
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where ^1/2, ^3/27 ••• ^ire the corresponding couphng constants. We adopted half-integer 
indices in order to foUow more easily the analogy with the microscopic formulas. The loop 



equations (|2.13| ) lead to the Virasoro constraints 



L^-ZN[t]=0 



with 



d d 



in > -11 



Ed 1 \ - , 



s — s' =n 



(5.2) 



(5.3) 



-s — s' =n 



Note that in the continuum limit the string equation is equivalent to the constraint defined 
by the operator L_i and the constraints ( |5.2| ) determine completely the partition function. 
The Hirota equations for the KdV hierarchy are formulated in terms of the vertex operators 



d 



s dt. 



(5.4) 



The loop equations and the bilinear functional equations in the continuum limit of the 
0(n) model are obtained by transforming the Virasoro generators and the vertex operators 
with the operator e^, where H is given by 



M- 



r(s + s' + i) 



d d 



2 frti « + ^' r(s + 1) r{s' + 1) dt, dt, 

s+s'>l 



(5.5) 



The parameter M measures the deviation from the critical point (we use the same notation 
as in Ijl^) so that 1/M is the inverse correlation length of the matter field. In the limit 
M 00 the deformation operator becomes the identity operator. 

Similarly, the statistical model associated with the adjacency matrix Gab is described 
by the deformation operator 



H 



M- 



r(s + s' + l) d d 

s' 



s + s' r(s + 1) r(s' + 1) at« at 



(5.6) 



acting on the set of coupling constants 



t = 



a = 1, r; s G + 



(5.7) 



The partition function is determined by the Hirota equations and r the subsidiary 
conditions (string equations) 



L^^.Z^[t] = (a = l,...,r). 



(5.8) 
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6. Concluding remarks 

Our construction represents an alternative nonperturbative formulation of 2D gravity 
suited for investigating the infrared properties of the theory. At very large distances the 
fluctuations of the matter flelds can be neglected and the theory decouples to r copies of 
pure gravity, where r is the number of points of the target space. At smaller distances the 
fluctuations become important, but still can be considered as a perturbation. A highly 
nontrivial fact, related to the speciflc realization of 2D gravity we are considering, is 
that the perturbation can be imposed by means of a differential operator acting on the 
coupling constants. As a consequence of this, the Virasoro constraints and the Hirota 
bilinear equations characterizing pure gravity survive in some form in the theory with 
matter flelds. 

The coupling constants ( |5.7| ) are related to the coupling constants associated with 
scaling operators as follows. The coupling constants are the coefficients in the expansion 
of the loop fields (always in half-integer powers of z) at z = 0, which is convergent in the 
circle \z\ < 2M. The coupling constants associated with scaling operators represent the 
coefficients in the expansion (in fractional powers of z) of the same loop fields at z — oo, 
which is convergent in the domain |2; + M| > M. 

Most likely the two realizations of the series, the one which we are considering and 
and the one as a matrix chain with interaction tvMaMa+i, have the same scaling limit. We 
therefore expect that both cases are described, in the continuum limit, by (r + l)-reduced 
KP hierarchy and the differential constraints generating a VFr+i algebra. 



The exact integrable structures associated with the Aj. models are discussed in ||T9[] . 
It is shown there that the canonical partition function of the Aq model is a r-function of 
the KdV hierarchy. Thus the KdV hierarchy appears in two cases: in pure gravity (the Ai 
model) and in the case of a gaussian field compactified at the Kosterlitz-Thouless radius 
(the 0{2) = Aq model). One can speculate that the partition functions of the Ar and 
Ar-i models are r-functions of the {r + l)-reduced KP hierarchy and that, more generally, 
the integrable structures associated with the ADE and ADE models are given by the 
exceptional hierarchies studied in [pO| . 
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